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Abstract: In this paper, we present a more detailed version of our previous work for three-particle 
correlations in quark and gluon jets [1|. We give theoretical results for this observable in the double 
logarithmic approximation and the modified leading logarithmic approximation. In both resummation 
schemes, we use the formalism of the generating functional and solve the evolution equations analytically 
from the steepest descent evaluation of the one-particle distribution. In addition, in this paper we include 
predictions beyond the limiting spectrum approximation and study this observable near the hump of the 
single inclusive distribution. We thus provide a further test of the local parton hadron duality (LPHD) 
and make predictions for the LHC. The computation of higher rank correlators is presented in the double 
logarithmic approximation and shown to be rather cumbersome. 
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1 Introduction 



The observation of quark and gluon jets has played a crucial role in establishing Quantum Chromody- 
namics (QCD) as the theory of strong interaction within the Standard Model of particle physics (2j[3). 
Jets, narrowly collimated bundles of hadrons, reflect configurations of quarks and gluons at short dis- 
tances EE). 

The evolution of gluon and quark initiated jets is dominated by soft gluon bremsstrahlung. Powerful 
schemes, like the Double Logarithmic Approximation (DLA) and the Modified Leading Logarithmic Ap- 
proximation (MLLA), which allow for the perturbative resummation of soft-collinear and hard-collinear 
gluons before the hadronization occurs, have been developed over the past thirty years (for a review 
see O). In the frame of high energy jets, one of the strikest predictions of perturbative QCD (pQCD), 
which follows as a consequence of Angular Ordering (AO) within the MLLA and the Local Partem 
Hadron Duality (LPHD) hypothesis [7], is the existence of the hump-backed shape (H of the inclu- 
sive energy distribution of hadrons, later confirmed by experiments at colliders like the LEP [|9j[T0l and 
the Tevatron ifTTl . Within the same formalism, the transverse momentum distribution, or k± -spectra of 
hadrons produced in pp collisions at center of mass energy yfs = 1 . 96 TeV at the Tevatron [12], was well 
described by MLLA [Q2J and next-to-MLLA (NMLLA) H41Q3J predictions inside the validity ranges 
provided by such schemes, both supported by the LPHD. Thus, the study and tests of enough inclusive 
observables like the inclusive energy distribution and the inclusive transverse momentum k± spectra of 
hadrons have shown that the perturbative stage of the process, which evolves from the hard scale or 
leading parton virtuality Q ~ E to the hadronization scale Qo, is dominant. In particular, these issues 
suggest that the hadronization stage of the QCD cascade do not affect pQCD predictions and therefore, 
that the LPHD hypothesis is successful while treating one-particle inclusive observables. 

The study of particle correlations in intrajet cascades, which are less inclusive observables, provide a 
refined test of the partonic dynamics and the LPHD. In lfl6l . the two-particle correlations inside quark 
and gluon jets were first computed at DLA. In H17II181 , this observable was computed for the first time at 
MLLA for such particles, whose energy or x (energy fraction of the jet carried away by one parton) stays 
close to the maximum of the one-particle distribution. In [19], the previous solutions were extended, at 
MLLA, to all possible values of x by exactly solving the QCD evolution equations. This observable was 
measured by the OPAL collaboration in the e + e _ annihilation at the Z° peak, that is for -^/i = 91.2 GeV 
at LEP Col . Though the agreement with predictions presented in |fl9l turned out to be rather good for the 
description of the data [20], a discrepancy still subsists pointing out a possible failure of the LPHD for 
less inclusive observables. However, these measurements were redone by the CDF collaboration in pp 
collisions at the Tevatron for mixed samples of quark and gluon jets [1 1 j. The agreement with predictions 
presented in |[T7l lT8 l turned out to be rather good, in particular for very soft particles (x <C 0.1) having 
very close energy fractions (xi « X2). A discrepancy between the OPAL and CDF analysis showed up 
and still stays unclear. That is why, the measurement of two-particle correlations at the LHC becomes 
crucial. 

By going one step beyond, in this paper we give predictions for three-particle correlations inside quark 
and gluon jets. This observable together with two-particle correlations can be measured in equal footing 
at the LHC. Such tests will provide further verifications of the LPHD for less inclusive observables 
and shed more light on the role of confinement in jet evolution. Further issues on the importance of 
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correlations versus single-particle distributions studies have been presented in H2 111221 . 
The paper is organized as follows. 

• in section[2]we recall the formalism of jet generating functionals and their evolution equations; 

• the kinematics and the process for the inclusive production of three particles inside the jet are 
specified in subsection 12 . 1 1 and 12.1.11 respectively ; 

• in subsection 12. 21 we obtain the MLLA exact system of integro-differential evolution equations for 
the three-particle correlations and in subsection 12.31 the single logarithms (SLs) contributions are 
obtained from the exact evolution equations; 

• in subsection 12.41 we obtain the DLA solution of the evolution equations and study the shape and 
overall normalization of this observable; 

• in subsection 12.51 the evolution equations are solved iteratively and the solution are expressed in 
terms of the logarithmic derivatives of the one-particle distribution and the two-particle correla- 



• in subsection 12.61 we finally give the analytical predictions which will be displayed in order to 
provide predictions for the Tevatron and the LHC; 

• in subsection 12.71 the hump approximation is applied to this observable; 

• in subsection 12.81 the region in x where the emission of three correlated particles becomes domi- 
nant is discussed; 

• in subsection 12.91 we give the analytical solution of the DLA four-particle correlator and show 
that including higher order corrections for differential higher rank correlators would become a 
cumbersome task; 

• in subsection |3j the predictions are displayed and the phenomenology is applied to the Tevatron 
and the LHC; 

• a conclusion summarizes this work; the appendices are written as complements of the main core 
of the paper. 

2 Formalism of the generating functional 

A generating functional Z(E, 0; {u}) can be constructed 11231 that describes the azimuth averaged parton 
content of a jet of energy E with a given opening half-angle 0; by virtue of the exact angular ordering 
(MLLA), it satisfies the following integro-differential evolution equation ||6l 



tions; 



dlnG 
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[z B {zp,Q;{u}) Z c ((l- z)p,Q;{u}) - Z A {p,Q;{u 



■») 



(1) 



2 



in (Q]), z and (1 — z) are the energy-momentum fractions carried away by the two offspring in the A — > 
BC parton decay described by the standard one loop splitting functions 11241 

*M{z) = C F ±±£ t ^{z) = C F l + { \- z) \ (2) 



$*®(z) = T R (z 2 + (1 - zf) , &W( Z ) = 2C A H—Z + + Z (1 _ z )j 



(3) 



C A = N C , C f = (N£-1)/2N c , T r = 1/2, (4) 

where N c is the number of colors; Z A in the integral in the r.h.s. of CO accounts for 1-loop virtual 
corrections, which exponentiate into Sudakov form factors. a s {q 2 ) is the running coupling constant of 
QCD 

, 0, 47T 

u s {q ) = — , (5) 

4iV c /3 ln-|— 
A QCD 

where Aqcd ~ a few hundred MeV's is the intrinsic scale of QCD, and 

is the first term in the perturbative expansion of the /3 function, n j the number of light quark flavors. 

If the radiated parton with 4-momentum k = (ko,k) is emitted with an angle with respect to the 
direction of the jet, one has (k± is the modulus of the transverse trivector k± orthogonal to the direction 
of the jet) k± ~ \k\@ w ko® w zEQ when z -C 1 or k± « (1 — z)E@ when 2 — >• 1, and a collinear 
cutoff fc^ > Qo is imposed. 

In (Q]) the symbol {u} denotes a set of probing functions u a (k) with the 4-momentum of a secondary 
parton of type a. The jet functional is normalized to the total jet production cross section such that 

Z A (p,e;« = l) = l; (7) 

for vanishingly small opening angle it reduces to the probing function of the single initial parton 

Z A (p,Q^0;{u})=u A (k=p). (8) 

To obtain exclusive n-particle distributions one takes n variational derivatives of Z A over u(ki) with 
appropriate particle momenta, i = 1 . . . n, and sets u = afterwards; inclusive distributions are gener- 
ated by taking variational derivatives around u = 1. We introduce the n-particle differential inclusive 
distribution, also known as parton densities, as J6] 

fnl $ n 

x\. . . x n D]i'(xi, . . . , x n ,Y) = Ei . . . E n ——— — — — Z A (kx,...,k n ,G;{u(k)}) . (9) 

6u(kx) ■ ■ -Ou(k n ) u=1 

Accordingly, we introduce the following notations for gluon and quark jets A = G, Q, Q 

41 (*) = y ••• y4 n) (7» + ln*), 4i = X! . ..^^(xi, . . . ,x n ,r), (10) 

which we will use hereafter; Xi corresponds to the Feynman energy fraction of the jet taken away by one 
particle In the case of three-particle correlations n = 3, the observable to be measured experimen- 
tally reads 

r (3) 41 
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2.1 Kinematics and variables 



The probability of soft gluon radiation off a color charge (moving in the z direction) has the polar angle 
dependence 

sinGde dsin(6/2) dQ 



2(1 - cos 9) sin(9/2) ' 
therefore, we choose the angular evolution parameter to be 

Qo sm(G/2) 
note that this choice accounts for finite angles 0(1) up to the full opening half-angle = it, at which 

Y e=7T = In— , 

where IE is the center-of-mass annihilation energy of the process e + e~ — > qq. For small angles (fTTT> 
reduces to 

r=!ta & e<<1 - £ = 3Ee- < 12) 

where Q = E@, defined as the virtuality of the jet, is the maximal transverse momentum of a partem 
inside the jet with opening half-angle 0. Moreover, we make use of variables known from previous 

works num, 

£ = \n—, y = \u^— — , A=- , (13) 

Xi &-QCD 

Li = In — , yj = In Xj ^ Q ° , rjij = In — , F = 4 + j/j + 77^. (14) 

Since ^ = dh fe 1 , y could also be used as the evolution-time variable in forthcoming quark and gluon 
jet evolution equations. Accordingly, the anomalous dimension, related to the coupling constant ([5]>, can 
be parametrized as follows 

7o V) = 2N C ^± => 7o V + y) = ] (15) 

such that, 

• for one particle 0, the denominator in ([T5l is simply £ + y + A, with ll26l £ = In |, y = In 

?7 = 0; 



• for two-particle correlation ||T9]|25], I + y + 7712, with I = In — , y = In X2 i? ei , 7712 = In a ; 

• for three-particle correlation, £ + y + 7713, with ^ = In — , y = In X3 /? ei , 7713 = 7712 + 7723 = In ^i. 

2.1.1 Integration bounds for three-particle evolution equations 

The production of three hadrons is displayed in Fig JI] after a quark or a gluon (A) jet of energy E, 
half opening angle ©o and virtuality Q = EQq has been produced in a high energy collision. The 
kinematical variable characterizing the process is given by the transverse momentum k± = zE®\ > Qo 
(or (1 — z)E@i > Qq) of the first splitting A — > BC. The parton C fragments into three offspring such 
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©1 



B (l-z)E 



Figure 1 : Three-particle yield and angular ordering inside a high energy jet. 



that three hadrons of energy fractions x\, x 2 and X3 can be triggered from the same cascade following 
the condition: 

e > 0! > e 2 > e 3! (16) 

which arises from the exact AO in MLLA [6]. In particular, the condition ©o > @i is kinematical rather 
than supported by the AO; it states that every collinear gluon is emitted inside the jet of half opening 
angle 0q. The two variables entering the evolution equations are z and ©i, such that 



xi<z<l 0<£<£l 
From (fT6l ) and the initial condition at threshold x 3 E@q > x 3 EQi > x 3 E@ 3 > Qq, one has 

Qo 



x 3 E 



< 01 < 0r 



< y < y3- 



(17) 



(18) 



2.2 From single inclusive distribution and two-particle correlation to three-particle cor- 
relation 

The evolution equations satisfied by (© are derived from the MLLA master equation for the generating 
functional ZA(u(ki)) ©. In this case, one takes the first second 5u{ ^ {k2) , and finally third 

Su(ki) Z Su(k 3 ) functional derivatives of ZA{u(ki)) over the probing functions u{ki) so as to obtain the 
system of evolution equations for 3-particle correlations. Following from dH), after applying the oper- 
ator su(ki) 3 6u{k 3 ) to k° tn members of the equation, according to (|9]l and (flOl together with the initial 
condition f7]), it is straightforward to get the coupled system of evolution equations 



r dz^ 9 q {z) [G^(z) + ( y Q^(l-z)-Q^+G^(z)Q 3 (l-z) + G 3 (z)Q^ 



+ G ( $(z)Q 2 (l -z) + G 2 (z)QY> + - z) + G x {z)Q 



)(2) 



(2), 



)(2) 
-23 



(19a) 



I dz—^ g (z) G®\z)-zG®+ Gg\z)G 3 (l -z) + G™ {z)G 2 {l - z) 



,(2), 



(2), 



+ G ( i{z)G l {\-z) +1 dz^n f m{z) (2Q( 3 \z)-GW)+2Q { $(z)Q 3 (1-z) 



j- 1 



a 

7T 



i( 2 )/ 



+ 2Q$(z)Q 2 (l -z) + 2Q^(z)Q l (l - z) . 
The l.h.s. of the equations d 19ab and dl9b|) can be written in the convenient form 



1(3) = A (3) _ Al A 2 A^ _ {A f) _ Ai a 2 )Az - {A™ - A 1 A Z )A 2 ~ (4 2 3 " A 2 A 3 )A 1 , 



(2) 



,(2) 



(19b) 



(20) 
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where A = G, Q, Q is the leading partem of the jet. Moreover, we have introduced the notations A\ 
4^(1), where 



(n) 



4 n i(l) =Xl...Xn£> (n) (xi,...,X n ,y), 



^l...n — ^l...nV 

for the sake of simplicity. The evolution equations for the single inclusive distribution and the two- 
particle correlation are written in ||T9l in the form 



Qy = fdZ^& q {z) 



G i, 



Q{1 - z)-Q) +G(z 



az — 

7T 



®>{z) G(z) -zG)+n f &(z) 2Q(z) - G 



(21a) 
(21b) 



and 



(Q(2) _ Ql Q 2 ) y = f\ z ^ <S>9 (z) Q m {z) + ( Q (2) ( 



+ (g 1 (z)-Q 1 )(q 2 (1-z)-Q 2 



1-z)- 
G 2 (z)-Q 2 )[Qi(l 



Q 



(22a) 



{G^-G 1 G 2 ) y = f dz^M(z) 

J xi n 



zG^) + (G l {z)-G l )[G 2 {l-z 
z) - Qi(z)Q 2 (z)) - (g®-G x G 2 



+ {2Q 1 (z)-G 1 ){2Q 2 (l-z)-G 2 



(22b) 



respectively. Making use of the equations (I21al2 lbb and (122al22bb . one can then construct the total 
derivatives [A^As] [(4? - A 1 A 2 )A 3 ] , [(4 3 - A]_A 3 )A 2 ] , [(4? - A 2 A 3 )A 1 ] as they ap- 

L -I y L J y L J y 

pear in (l20l >. which are to be subtracted, term by term from the system of equations dl9a|19bj ). Therefore, 
we get the equivalent system for the three-particle correlations inside quark and gluon jets: 

Q y 3) = I dz^$9{z)[GW{z) + (Q( 3 Xl-z)-QW) (23a) 



j-i 



+ (Q? 2 \l 



Q { 3) (G 3 (z) - g 3 ) + (d§{z) - q§) (Q 3 (i -z)-q s 



+ (gg>(l - 2) - Qg) (G 2 (z) - Q 2 ) + (g$(z) - Q<?) (Q 2 (l - z) 



0.2 



Q23 ) (Gi(^) - Qi) + (G 



.(2), 
T 23 I 



Q^KQiCl-zJ-Qi 



+ ((Qi - Gi(*)) (Q 2 (l - z) - Qa) + (Q2 - G 2 (z)) (Q a (l - z) - Q0) Q 3 
+ ((Qi - G x {z)) (Q 3 (l - 2) - Q 3 ) + (Qs - G 3 (z)) (Qa(l - z) - Qj)) Q 2 
+ ((Q2 - G 2 (z)) (Q 3 (l - z) - Q 3 ) + (Q 3 - G 3 (z)) (Q 2 (l - z) - Q 2 )) Q x ] , 



Gg ) )(G 3 (l- 2 )-G 3 ) 



+ (Gg(z) - Gg>) (G 2 (l - z) - G 2 ) + (<?§?(*) - Gg } ) (Gi(l - z) - d) 
+ (Gi - Gi(z))(G 2 (l -z)- G 2 )G 3 + (Gi - Gi(z))(G 3 (l - z) - G 3 )G 2 

.)-G( 3 )) 



+ (G 2 -G 2 (z))(G 3 (l-z)-G 3 )G 1 ]+ / tfe-^n/*«( 



+ 2(QS 2 2 ) (z)-Gg ) )(g 3 (l 
+ 2 (q ( 3(z)-G?/)(Q 2 (1 



G 3 ) + (2Q 1 (z)Q 2 (z)-G 1 G 2 )G 3 
G 2 ) + (2Q 1 (z)Q 3 (z)-G 1 G 3 )G 2 



(23b) 
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+ 2 (q$(z) - Gf 3 } ) (Qx(l - z) - Gi) + (2Q 2 (z)Q 3 (z) - G 2 G 3 )Gi 

+ (Gi - 2Qi(z))(2Q 2 (l - z) - G 2 )G 3 + (d - 2Qi(z))(2Q 3 (l - z) - G 3 )G 2 

+ (G 2 - 2Q 2 (z))(2Q 3 (l - z) - G 3 )Gi] . 

The system of evolution equations (I23al23bb . which appears as a consequence of the exact AO in intra-jet 
cascades, provides the complete theoretical picture of the three-particle correlations as a function of %i 
and the characteristic hardness of the jet Q; this is the first new result of this paper. However, since these 
equations could only be solved numerically, we will extract the SLs contributions O(y^aJ) in order to 
provide an approximated analytical solution in the following. 



Qf= / dz^& q (x)G®(z), (24) 
Gf = fdz^(l - z)*>{z)GP>{x) + f dz^n f $l(z) [{2Q^ - G< 3 >) + 2 (flg - G§>) (25) 

J X~\ <J X 1 



2.3 Approximate evolution equations 

Let us start with equation (I23ab . We proceed to cast all SLs contributions corresponding to hard-collinear 
parton splittings in the shower. In the hard parton fragmentation region one has z ~ (1— z) ~ 1, such that 
the second contribution in (|23ab can be approximated through a Taylor series for In z ~ ln(l — z) <C 
written in the appendix |A] Therefore, one obtains the simplified system of evolution equations 



x (Q 3 - G 3 ) + (2QiQ 2 - GiG 2 )G 3 + 2 - G\1> ) (Q 2 - G 2 ) + (2QiQ 3 - GiG 3 )G 2 

+ 2 (Q® - Gg } ) (Qi - G x ) + (2Q 2 Q 3 - G 2 G 3 )Gi + (Gi - 2Qi)(2Q 2 - G 2 )G 3 
+ (Gi - 2Q0(2Q 3 - G 3 )G 2 + (G 2 - 2Q 2 )(2Q 3 - G 3 )Gi] , 

where we have kept all terms of order 0(^/a~ s ), which contribute to MLLA. In addition, from the DLA 
relation Z A = Z^ a/Nc (271, and Eqs.dSED, one has the useful expressions for the single inclusive 
distribution, two- and three-particle correlations: 

^ ] - % g,3) + l(f - Ok* + + c - G + % (f - 0(1 - 2 

x GiG 2 G 3 , ' (27) 



which in turn can be replaced in (I25l l. The two expressions written in (1261 ) are known from previous 
works at DLA H161I27II . while (|27T ) will be used for the first time in this context. After integrating over 
the regular part of the splitting functions (|2), ((3) and ®, one obtains the integro-differential system of 
equations (7/13 = rjxs + ^23), 

Qf = ^fJ^l 2 o(i + y3)G^(£,y 3 ;m3) - ^7oVi + l/3)G< 3 >(4, y 3 ; r? 13 ), (28) 

Gf = f 1 d£^(£ + y 3 )G ( - 3 \£,y 3 ;nx 3 ) - a 1 i(£x+y 3 )G ( - 3 \£x,y3;rix3) + (« - &bg(4 +lfe)(29) 

JO 

x (G^(4,y 3 + r/ 23 ;r/i 2 ) - Gl(£i, 2/3 + 7713) G 2 (4l + ^12, y3 + ^23)) G 3 (^i +7713,7/3) 
+ (4, 2/3; 7713) - Gi(^i,t/ 3 + r/i 3 )G 3 (^i + 7713,2/3)) G 2 (4 + ??i 2 , y 3 + r/ 23 ) 
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+ {p 23 \h + m2,y3;m3) - G 2 (i x + m2,y3 + m3)G 3 (£i + 7713,2/3)) Gi(h,y 3 + rj 13 

+ {a-c)j%(e 1 + y 3 )G l (£ 1 ,y3 + ms )G 2 (h + V12 , V3 + V23 ) G 3 (l x + vi3,ys), 



with the following hard constants, 

1 



a{n f ) 
b(n f ) 

c(n/) 



4N C 
1 



4/V r 



11 4 C F 
—N c + -n f T R 1 - 2— 
3 3 \ N r 



n/=3 



0.935, 



11 



-iV, 



;n f T R 1-2 



iV, 



rtf =3 



0.915, 



HiV c + t nfTR (1-2Q 

3 3 1 V iVr 



71/ =3 



0.917, 



(30) 
(31) 
(32) 



where nj = 3 corresponds to the number of light active flavors of quarks u, d, s. As an example of such 
procedure, one could write the example, 



a(n f ) 



The first integral terms of the equations in (I281 l and (1291 are of classical origin and therefore, universal. 
Corrections oc — |, a, (a — b) and (a — c), which are 0(- s /al) suppressed, better account for energy 
conservation at each vertex of the splitting process, as compared with the DLA. Notice that the form of 
the quark initiated jet equation (|28T ) is universal at MLLA (see (l80l ) and (l82~1in the appendix lA.ll for the 
single inclusive distribution and two-particle correlation respectively), that is, invariant with respect to 
the number of particles considered in the cascade. In the equation for the gluon initiated jet (1291 ), the 
first and second constants a(rif) and b(nf) were obtained in the frame of the single inclusive distribution 
and two-particle correlations respectively ITT71IT81 . The third constant c(nj) appears in this paper for 
the first time for the three-particle correlation. In particular, notice that a certain recurrency shows up 
in the coefficients combining the colour factors (— l) n_1 ^1 — 2^^ , as a function of the number n of 
particles considered in the shower. 



2.4 DLA solution of the evolution equations 



In this subsection we compute the leading order DLA contributions in order to provide general features 
concerning the the shape and overall normalization of three -particle correlations. This procedure is 
equivalent to cast the leading order (LO) solution of the equations (1281291) . We differentiate (l28l ) and 
( [291) with respect to such that after setting hard corrections oc 3/4, a, b, c = 0, the MLLA evolution 
equations are reduced to the new DLA compact differential equation 



i(3) 



iy N c 



Ca„2 g(3 ) 



■7o< 



(33) 



with 



i(3) 



c 



(3) 



c 



(2) 



,(2) 



c 



(2) 



AiA 2 A 3 



ty 



(34) 



after having set = Cj? AiA 2 A 3 for the three-particle correlator and A^ 



(2) 



C^ AiAj for the two- 



particle correlator. We fix the anomalous dimension to the characteristic hardness of the jet Q « EQq 
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(7o(^©o) = const) and solve this equation iteratively by derivating the r.h.s. of (|34l > with respect to 
and y, such that the solution of (l33l) reads 



4 3 L-1)-(C£-1)-(C£-1)-(C£-1] (35) 



C A 2 + A 12 + Aw + A 23 Cl 2 + A 12 + Ai 3 + A 23 ' 

which have been written in terms of the logarithmic derivatives of the one-particle spectrum, 

Ay = 7cT 2 [$A i fil>A i ,v + ^A^A^l) , jlAiJt = j:^f> ^y = (36) 
and the DLA two-particle correlator EQiD (for a review see also [28 ]) 

The dot over differentiates the DLA correlators from the MLLA correlators obtained below. In DLA 
however, since the single inclusive distribution satisfies Q = j^G IT271 . one has 

That is why, we will use the much simplest notation ipc z ,i = "<pGi,y = ipi,y It is worth giving the 
order of magnitude of some quantities that will be considered in forthcoming calculations. In DLA, the 
one-particle inclusive distribution can be written as Ai(£,y) oc exp (2joy/£y) asymptotically for fixed 
running coupling 70 = const ll27l . Though the solution with fixed coupling constant provides general 
features of the single inclusive distribution, it is not enough for the description of a more realistic picture 
at colliders. However, from its simplicity, it can be used to give the order of magnitude of terms involved 
in the solution of the DLA and MLLA evolution equations. Therefore, making use of d36l ), one has 

^A h i = O(7o), ^ Ai ,y = O(7o), = °(7o)> ^A^y = C(7o), = C(7q), (38) 

\j = 0(1), A tj/ = 0( 7o 2 ), K ihy = 0( 7o 2 ), (39) 



where V'A 4 ,«» ^Uj,% an d ^Ai,yy are double derivatives of 1/}^ = In Ai(£,y). The DLA solution ( T35l) 
describes the following picture: the first term (= —1) in the l.h.s. translates the independent or decor- 
related emission of three hadrons in the shower like depicted by Figj2^. After inserting the two-particle 
correlator (l3"7T i in the l.h.s. of (|35T >. terms oc ^ correpond to the case where two partons are correlated 
inside the same subjet, while the other one is emitted independently from the rest like in Figf2j). Next, re- 



placing (I37T ) in the r.h.s. of (l35l l one obtains a contribution oc described by Figf2£, where two partons 



C 



A 



N 2 

are emitted independently inside the same subjet. The last term oc depicted by Figj2ti, involves three 

G A 



particles strongly correlated inside the same partonic shower and corresponds to the cumulant of genuine 
correlations. Actually, this interpretation has been given after computing the color factors of such Feyn- 
man diagrams describing the process, normalized by C\ in the end. Notice that diagrams displayed in 
Figf2fc and Figj2ji present the same color factors but different Lorentz structure. In both cases, the DLA 
strong AO 6>6'> 6" and strong energy ordering x\ 3> x 2 3> x 3 are necessary conditions satisfied 
by (291. 
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A 




Figure 2: Three particles emitted inside the shower with color factors for the square of the amplitudes: 
C\, C\N C , CaN£ and CaN% for a, b, c and d respectively. 



Performing the steepest descent evaluation of the DLA single inclusive distribution from an integral 
representation, which was written in Mellin space in the form II161I27I . 



G(£,y) = (£ + y + X) 



dojdv 
(2^)2 



ds ( ijj(v + s) 



o v + s \(oj + s)v 



1/00 (w-«/) 



e~ Xs , Q = ^fG. (40) 



and which accounts for the running of the coupling a s , the energy of most particles inside the jet was 
proved to be close to the maximum of the distribution, which shapes like a Gaussian in this region ll27l . 



Ai(£i,Y) ~ exp 



Y 



j ^max 



(41) 



% y 3 /2 

From this method [16], the expressions of the logarithmic derivative of the one particle distribution were 
written as, 

A,l(tH, n) = 7oe M % ipi, y ((i, v) = 7 e~ M \ (42) 
such that Ajj and the correlator were given in the form ||T6l , 

N r 1 



' (2) 

Aij = 2cosh(fii — Uj)j C y A } _ = 1 + 



(43) 



Ca 1 + 2 cosh(/ij — fij) 

respectively, where (/Xj, were related to (4, Ui) through the 2x2 non-linear system of equations |[T6l . 
Ui — ii (sinh 2/ij — 2^j) — (sinh 2fj — 2i/j) sinh^j sinh/Uj 



Vi+£i 2(sinh 2 /ij - sinh 2 ^) \f\ V^i + Vi + A 

Therefore, the DLA three-particle correlator reads in this approximation 



i(3) 



5(2) 
Mia 



+ 



+ 



Nr. 



1 + c 



;(2) 

'Al3 



1 + c 



i(2) 
M23 



1 



(44) 



(45) 



(2) 



J(2) 
Ml 3 



;(2) 



2Ca 1 + cosh(/xi - /i 2 ) + cosh(/ii - /z 3 ) + cosh(^ 2 - ^3) 



Ni 



1 



2Cj| 1 + cosh(^i - /j, 2 ) + cosh(/ii - // 3 ) + cosh(^ 2 - ^3) ' 



with extracted from (|43T >. Taking 
approximation (see appendix IC.2I ) 

Ay « 2 + 9 



|<g; cr oc y 3 / 2 for i = 1,2,3, one has in this 



Y 



2 + 9 



ln(xj/x 4 ; 



i 2 



ln(Q/Qo) 



(46) 
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so that, 



A12 + A13 + A 



^23 



6 + 9 
6 + 9 



Y 

ln(x 2 /xi) 
HQ/Qo) 



+ 9 

2 

+ 9 



Y 

\n(x 3 /x 1 ) 
HQ/Qo) 



+ 9 

2 



Y 



+ 9 



m(x 3 /x 2 ) 
HQ/Qo) 



(47) 



Therefore, the shape of the three-particle correlator can be expected to be quadratic as a function of the 
difference — £j), as for the two-particle correlator. Thus, the correlator is strongest when particles 
have the same energy x\ = x 2 = x%. 

Moreover, the decreasing behavior of the correlator as one parton gets much harder than the others 
xi 3> xj shows that QCD coherence effects dominate this region of the phase space as interferences 
between such gluons occur. New kinds of contributions like the one in the first term of the r.h.s. of (l35l l 
appear in this context. 

The overall normalization of the n-particle correlator is fixed by that of the same rank multiplicity- 
correlator determining the multiplicity fluctuations inside the jet lfl6ll . 

(n(n - 1) . . . (n- k+ 1)) 



C^\x±, ...,Xk)< 



(n) 



Then, one has 



C^ 2) (x 1 ,x 2 ) 



1 < 



Nr. 



Cy(x 1 ,x 2 ,x 3 ) - 1 < — + 



(48) 



3<V ~ A 

These bounds can also be obtained by setting A(xi, Xj) = 2 (for x% = Xj) in (|37T > and (l35l) respectively. 
Since DLA neglects the energy balance, it is not realistic and does not provide the real physical picture 
of any jet process in the frame of jet calculus. 



2.5 Iterative solution of the evolution equations 

As we can see, the computation of three-particle correlations requires a mastering knowledge of the 
one-particle inclusive energy distribution and two-particle correlations. The behavior of the two-particle 
correlators as shown by these solutions was proved to be quadratic as a function of —lj) and increasing 
as a function of (£i + £j) like in the Fong-Webber approximation H17II18II . However, the solutions (192193b 
(see appendix IA.1I ) were shown to better account for soft gluon coherence effects, by describing the 
flatting of the slopes as + lj) increases. In ||25l , the solution was obtained by the steepest descent 
evaluation of the spectrum Gi(£, y ), while in [19 ], the evaluation was performed by taking the expression 
of Gi(£,y) given by d89l in the appendix IA.1I In [19], the solution of the evolution equations for two- 
particle correlation were obtained from the differential version of the equations (1901911 ) over £ and y 
written in the appendix IA.1I Therefore, in this subsection, we will make some transformations in order 
to simplify this cumbersome task without adding further information. In the appendix IA.1I we briefly 
summarize what should be known in order to complete the solution of the evolution equations for the 
three-particle correlations. 

Differentiating d28l and d29l with respect to "f , one has the differential system of evolution equations 



11 



for three-particle correlations, 



G (3) 



N r 4 iV r V 



&$ = ^ -a 7o 2 (Gf > -/3 7 2 G (3) ) + (a - &) 7o 2 { [( 



g! 2 2 } - dG 2 ) G 3 



(49) 
(50) 



+ (Gg- g ,xv :: j (v 2 [a 

,(2) 
r 13 



+ G 



G1G3 I G2 + I G 



,(2) 
T 23 

,(2) 
T 23 



G2G3J Gi 
G2G3 Gi 



-A)75 



G 



(2) 
12 



G1G2 G3 



} + (a-c) 7o 2 [(G 1 G2G 3 )£-/3 7oG'iG2G3] 



which is written in this paper for the first time. The equation d50l is self-contained and can be solved 
iteratively like £S]>. For this purpose, one sets = C§ ) 123 G 1 G 2 G 3 and gJ 2) = C^dGj in the 
left and right hand sides of (l50l . such that the solution obtained in the appendix |B] can be written in the 
compact form 



C. 



(3) 
G12 



1 



^(2) 
'G12 



:)^ 2) + U' 



-(2) 



1 W?( 2 ) + F (3) 
1 ' r 23 ^^123' 



where, 



(2) 



AT, 



1 + 



(2) 



(2) 
G 



(3) 



N. 



(3) 
G 



123 



(3) 
G 



with 



D 

N. 



(2) 

(2) 
G 
(3) 
G 



1 - b (Vi/ + V>2,£ + 4>3,e - A)7o) " o& + (a- b) X \ j + d J + $ 

2 + A12 + A13 + A 23 + a& + 2a/3 7o + d + e 2 , 



ei - £2, 



(51) 



(52) 



(53a) 
(53b) 



(VV + + ^ - /3o7o) - «C£ + (a - b)( X ? + xf + X? ) + (£f 2 + ^ 12 )(53c) 



(3) 
G 



+ (£l 13 + 4 3 ) + (£l 23 + <f ! ) 



ei - £2, 



(2) 
G 



2 + A12 + A 13 + A 23 + a& + 2a/3 7o + £ i + e 2- 



(53d) 



The solution (f5TT > can be checked to recover the DLA result (l35l inside a gluon jet, that is for Ca = iV c . 
Since DLA neglects recoil effects at each splitting inside the cascade, one should expect the DLA three- 
particle correlation to be much larger than MLLA predictions and therefore to overestimate the data. We 
introduce the following notations and give the order of magnitude of each contribution following from 
(EH) and 



C = lnC 



(3) 

G123 ' 



c. 



(3) 

Gl23,<? 



c, 



(3) 
G123 



O(7o), Cy 



C 



(3) 

Gi23,y 



c 



(3) 



O(7o 2 ), 



l 3 



c. 



4(2) 



c (2) 

G»j ,i 

a{2) ■ -,i 



5(2) 



O(7o 2 ), 



(1 



(1 

e 2 



7o 2 ^ 



XI {4>l,y + 1p2,y + fo,y) + Xy {^1,1 + 4>2,t + i>3,t) = O(7o) 

O(7o 2 ), 



^2 [XtXy + 
7o 



-2 (<Xy + Q,y) 

7o 



O(7o 2 )- 



(54a) 

(54b) 

(54c) 
(54d) 
(54e) 
(54f) 
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The solution of the gluon evolution equation for the correlator can be either obtained numerically by 
solving (l50l > or by performing the evaluation from the previous solution (I5TI ). However, in this paper, 
we will directly compute the solution (f5TT > from the steepest descent method introduced in ll25l and 
make some approximations in subsection 12.61 Accordingly, the solution of d49l ) is also obtained in 
the appendix E by setting = C { q\ 23 QiQ 2 Qz and Q^f = C^QiQj in the l.h.s. of @) and 
= Cq , G1G2G3 in the r.h.s. of the same equation, such that, 



where, 



with 



C. 



(3) 
Q123 



c, 



(2) 



,(2) 
'Ql3 



l)^ ( 3 2) +(^ 



,(2) 



123' 



F. (2) 



N, 



1 + 



(2) 

Q % j 



D 



(2) ' 

Q 



F 



(3) 
123 



N. 



(3) 



(3) ' 



(55) 



(56) 



Qij 

,(2) 



iV. 



(3) 

Q 



: gi + ^_ gl _g 2j 



L>g } = A 12 + A13 + A 23 + -%f + ei + e 2 , 



^ C (3) 



+ V>2,£ + + - /3o7o) 



GIG2G3 
Q1Q2Q3 



+ (£ 2 + 3 



12\ 



+ (£f + ^) + (£f + ^ d )-e 1 -e 2 , 



D 



(3) 



g)=A 12 + A 13 + A 23 + E% + ^+^ 

• 7oVi 



(57a) 
(57b) 

(57c) 
(57d) 



where one find the list of corrections, 



C = lnC 



(3) 



C. 



(3) 

Q123, 



c, 



(3) 
Q123 



0(7o), C« 



C. 



(3) 

?i23,y 



C. 



(3) 
Q123 



O(7o 2 ), 



c. 



(2) 



c. 



(2) 



c, 



4(2) 



Qjj,y 

j(2) 

Qij 



Oill) 



ei 
e 2 



O(7o 2 ), 



7o 2 ^ 

7o v 
l r 



7o 2 ^ 
1 



Ce&Quv + ^Q 2 ,y + $Q 3 ,v) + CvWw + + ^W) = C(7o), 



O(7o) 



-2 ( CeCw + 
7o 



o(7o 2 ; 



(58a) 

(58b) 

(58c) 
(58d) 
(58e) 
(58f) 



The order of magnitude of these terms follows from (1381 and (|39b . Setting all corrections to zero, one 
recovers the DLA solution (|35T > for Ca = Cf- The solutions dFTT ) and d55l ) of the evolution equations en- 
tangle corrections of order 0(70) and 0(7q), which are MLLA and NMLLA respectively. Furthermore, 
every term in (1511 and (|55T ) can be associated to a Feynman diagram of Figf2]as was explained in subsec- 
tion !2.4l The functions F^ 3 and F^l in (IBTT) and d55l) correspond respectively to the cumulant of genuine 
correlations associated to the process displayed in FigfJJand Fig|2ji. These contributions, d54aH54fb and 
d58aH58f| > are small corrections arising from the iterative solution of the evolution equations because one 
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takes the derivatives over the functions C = In , C = In C>P and y jj = In , y lJ = In for 
both quark and gluon jets. For the evaluation of such corrections one needs to take the DLA expressions 
of Ca 123 and . written in (|35T ) and (I37T ) respectively. 

2.6 MLLA approximation and evaluation by the steepest descent method 



In 1119 1, the exact solutions of the two-particle evolution equations were compared with the MLLA so- 
lutions from the steepest descent method for the one particle distribution. The agreement between both 
approaches was successful and made possible the fast computation of the correlators from the steepest 
descent. That is the reason for in this paper, we limit ourselves to this method. Making use of the ratio 
d87l ), it is easy to demonstrate that, 

V>Q/ = V>£ + 0(7o)> fe = ^ + 0(7o), A^ a A ij + 0( 7 o 2 ). (59) 

Dropping corrections of order 0(7q), which go beyond the MLLA approximation, we obtain for the 
gluon jet 

p (2) mUa l-b Q M + 4> 2 ,e + *03,l) + Cl ~ ^1 (60) 

ij 2 + A12 + A13 + A 23 + ei 



p (3) mUa 1 ~ C + ^ + fa, t ) + g| 2 + gf 3 + gf - Cj 

2 + A 12 + A 1 3 + A 23 + 6 1 (61) 



and for the quark jet 



^,(2) m/Za j _^ £1 gl ,g 2 ) 

« 3 + A12 + A13 + A 23 - a + W + foj) + h ' 



p{ 3) mUa N 2 C Cg^ [1 ~ a (Vy + V> 2 ,€ + jMj + fg + gf 3 + fg 

123 ' C\ 3 + A i2 + A13 + A 23 - a (Vv + i>U + + ei 



(63) 



The subtracted terms oc —a in the denominators of (l62l and (l63l appear after having replaced (I87T ) and 
(l88l ) in (I57bl ) and (I57cb respectively. Such simplified expressions are useful for the steepest descent eval- 
uation that proved successful while describing the single inclusive distribution and two-particle correla- 
tions in EBTl . Except the MLLA corrections e\ and , all the other corrections and functions appearing 
in the solutions of the evolution equations were obtained in [25 ], which will allow for the straightforward 
computation of the three-particle correlators in quark and gluon jets. We write some of these formulas 
for the evaluation in the appendix [C] Integrating the equation ( T8TT > over u y", the solution for the single 
inclusive distribution is given by the following integral representation in Mellin space ESI , 



G{1 , y ) = {t + y + X) J J j, ^ J ^ e ->. (64) 

The integral representation (l64l was estimated by the steepest descent method at small x -C 1 and 
high energy scale Q » 1; the approached solution was compared with the exact solution (|89l (see 
the appendix IA. lb in the limiting spectrum (A = 0) and beyond (A / 0). In particular, (l64l was also 
demonstrated to be equivalent to (l89l for A = lfT9l . The agreement between the approached and 
exact solutions turned out to be good, such that the following expressions of the approached logarithmic 
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derivatives from the steepest descent method were suited for the evaluation of the two-particle correlators 



1 



1 2 



e w <5(/Ltj, Vi) — tanh V{ — tanh Ui coth /ij I 1 + e w Q(/Ltj, ^ 



(65) 



7 e 
1 



1 2 

2° 7 ° 



/3o7o 1 + tanh z/i 1 + if )+C{jJn,Vi) [l + ePQfavM + 0( 7 g), 



2 + e ^Q^ii Vi) + tanh — tanh i/j coth /ijf 1 + e ^Q{^i,Vi) 



/3o7o 1 + tanh^ l + if^,^) -C(/ii,^) l + e-^Q(/ii,^) + 0(7o), (66) 



where the functions Q(fii, Vi), C(fii, v^) and if (/ij, i/j) are defined in the appendix ICl The term oc a in 
(|65T > and (|66T > accounts for energy conservation while that oc /?o accounts for the running of the coupling 
a s . The variables (m, Ui) are related to (£i,yi) through the same 2x2 non-linear system of equations (04). 
After inverting (|44T i numerically, Hi(£i,yi) and Ui{ti,Ui) can be plugged into ( [651 and (|66l ) so as to get the 
logarithmic derivatives of the single inclusive spectrum as a function of the original kinematical variables 
ti and yi as it was done in EBll . The MLLA two-particle correlators involved in (IBTT ) and (|55T ) are (11081 ) 
and d 109b and are written in the appendix |C] These expressions have been taken from reference ||25l . 

Corrections C\ , Ci an d ei,ei are new for three-particle correlations. Such expressions are explicitly 
written in the appendix IC.ll from the steepest descent evaluation of the single inclusive distribution 
(l64l ). They are small and decrease the three-particle correlator for 4 ^ £j, that is when one parton 
is much harder than the other. Notice that the steepest descent method constitutes the only way for the 
disentanglement between MLLA O{y/ol) and NMLLA 0(a s ) corrections appearing in the solution of 
the evolution equations for the two and three-particle correlations. It makes also possible to distinguish 
between corrections following from the energy balance and the running effects of the coupling constant 
a s . Finally, this method also allows for the application of the hump approximation or Fong-Webber 
expansion of the solutions with MLLA 0{y/al) accuracy H171I18I . 

In this frame, the role of MLLA corrections should be expected to be larger than for the two-particle 
correlations. Indeed, higher order corrections increase with the rank of the correlator, which is known 
from the Koba-Nielsen-Olesen (KNO) problem for intra-jet multiplicity fluctuations H28U30LI3TTI . For the 
2-particle for instance one has oc —b(ipi t £ + ^2,e) an d for the three-particle correlator one gets the larger 
correction oc -c(V>ij + Vv + V>3/)- 



2.7 Hump approximation 

From the steepest descent evaluation introduced in [25], near the hump of the single inclusive distribution 
| i — Y/2 |<S g oc Y 3 / 2 for i = 1,2, 3, corrections £f , ^ and e±, e\ could be written in the symbolic 
form (see appendix |C.2| ), 



e^e^(^) 70 + O(7o z ), (67) 



2 / n b \ 2 / „ « \ 2 



■ (^) 70+ fc^J 70+ (^j + (68, 



such that both can be neglected ^ rj 0, e\ « in this approximation, like tf£ was also in EBl . In the 
appendix IC.21 following from the steepest descent method, the expressions of (|53aH53d| ) are given and 
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d57aH57d| ) expanded in ,</a~ s . In particular, the expressions (|128et and ( 128g| l, after being expanded in 
70, can be demonstrated to recover the Fong -Webber results for the two-particle correlations lli"7lli~8~l . 
Replacing the expressions (I128al 128j I into (15 1 152b and (155 15 6b . one finds those for the three-particle 
correlators in the Fong-Webber approximation ITT71H81 . This solution will be compared with that from 
d6T| ) and ( f63b after making use of d65l and ( f66b in subsection [3] 



2.8 From two to three-particle correlations in the small x region 

In |fi"9ll , the sign of the two-particle correlator (C^ — 1 > 0) was studied as a function of x in the region 
of the phase space where the two partons (hadrons after assuming the LPHD) are strongly correlated. 
From the previous inequality, it turned out that two partons with £i > 2.6 (xj < 0.07) at LHC energy 
scales (i.e. Q = 450 GeV, see subsection [3]> are correlated as they are emitted from the same cascade 
following the QCD AO. Asymptotically Y -> oo, one has £( > 4.5 (x; < 0.011). 

(3) 

For three-particle correlations we study the sign of the cumulant of the genuine correlator F 12 g > and 
determine the approximate region in x where diagrams displayed in FigfT|and Figf2]l become dominant. 
One has, 

1 - C (W + fa* + M + t\ 2 + £l 3 + if - €1 > 0. 

However, corrections £j J , ei have been shown to be negligible and to vanish for particles having the same 
energy momentum. Thus, we rather study the sign of 

1 - C (lpl £ + 1p2t + lp 3 j) > 0. 



Making use of ipg = <jqJ f = jq\/ ±j± for the sake of simplicity, one has, 



Y — I „ M 9c 2 

W 
Y 



l-3c70A/^>0^^ >T ^, M = — = 10.1. 



Thus, for LHC energy Y = 7.5, the value of l{x) where the cumulant becomes positive turns out to 
be t > 4.3, which in x corresponds to x < 0.014. Asymptotically Y — > oo, one has li > 10.1 
(Xi < 4.1 x 10~ 5 ). Therefore, there exists a range in x where the observable C^ 3 is dominated by the 
emission of two correlated partons emitted independently from the third one, that is 0.014 < x < 0.07 
for diagrams Figf2j) and Figfl};; for x < 0.014, the process will be dominated by three particles emitted 
from the same partonic cascade following the QCD AO described in Figf2jl. Asymptotically Y — > oo, 
one has 4.1 x 10" 5 < x < 0.011 for diagrams FigfJh and Fig©, and x < 4.1 x 10~ 5 for Figj2jl. These 
values will indeed justify our choices for the representation of the three-particle correlations as function 
of (a?i , X2, Xs) in subsection [3] 



2.9 Beyond three-particle correlations 

It is worth reminding that the LPHD hypothesis has also been confronted to multi-particle factorial 
moments up to the 5th order in the experimental studies of ep and e + e~ collisions at HERA f34l and LEP 
|[35l respectively, where it was found that the LPHD hypothesis faces difficulties when it is applied to soft 
multi-particle fluctuations. In this work the studies are carried out by using the momentum and transverse 
momentum cuts in order to test the MLLA soft limit calculations [33]. The theoretical computation of 
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multiplicity correlators or multiplicity fluctuations [n[n — 1) . . . (n — k + 1)) was performed in 11321 at 
MLLA up to the rank k = 5 of the correlator. 

However, performing these calculations for higher rank differential inclusive correlators, related to the 
previous ones by the integral 

becomes rather cumbersome. As an example, in this subsection, we display the DLA equation and 
solution of the 4-particle correlator. The DLA equation reads, 

^1234 — ^70^1234' <- oy J 

where A has been defined in the appendix [D]in (11291) . The solution of (l69l with the definition of A (11291 ) 
reads, 

C f.^ Hl (d«) + | H> (^) + g , +Au+Aa g ;^ + A M+ A M - <™> 

where the functions i?2 and i?3 are written in the appendix |D]in (1 1 30b . (1131b and (11321 ) respectively. 
The solution (1701 ) can also be interpreted in terms of Feynman diagrams contributing to the emission of 
four hadrons inside the jet. Accordingly, the term oc ^ correspond to the case A — > 12(34) where 
two offspring are correlated while the other two are emitted independently; as a consequence it depends 

TV 2 

only on the two-particle correlator. The second term oc 7^- is associated to the cases A — > (12) (34) 
and A — > (123)4, which translates into either emitting two sub-jets with two-particles correlated within 
each, or emitting three correlated partons like in FigJUwith another independent emission. Finally, the 
term oc after setting #3 = 1 + . . . corresponds to the full correlated emission of four offspring inside 
the same shower. The inclusion of SLs corrections to (1701 would be cumbersome and stays beyond the 
scope of this paper. On the other hand, the computation of differential higher order rank ( k) correlators 
at MLLA would imply the failure of the perturbative approach because of the increasing size of higher 
order corrections oc (ipi t e + ■ ■ ■ ipk,e) = Hence, for higher order k correlators, the small x range 

where MLLA predictions stay valid gets reduced even at high energy scales, such that (see subsection 

M k = — * 4 > 



with 

Cfc 



4iV, 



Po ' 1 + ^ 



£«r. + (-l).Wl- 2 f 



3 Predictions for the LHC and phenomenological consequences 

In this section, we perform theoretical predictions for three-particle correlations for the LHC. We display 
the MLLA solutions (I5TI ) and (l55l) of the evolution equations (l50l and d49l respectively. We compare 
the DLA solution of the evolutions equations from section |Z41 with the MLLA solution from the steepest 
descent evaluation of the one-particle distribution in subsection 12.61 and the solution from the hump 
approximation in 12.71 Thus, 
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• the DLA solution is computed by plugging d43l into d45l ); 



the MLLA solution from the steepest descent will be displayed by substituting the MLLA two- 
particle correlators (11081 . (1109) and the functions (|60l . (l6il . (l62l and d63j into d2B and ^55]> for 
gluon and quark jets respectively; 



• the MLLA hump approximation will be displayed by plugging (1128al) -( 128j ) into (1521 and d56l ) 
and finally (TST]) and d55l ). 

In particular, the computation of the DLA and MLLA solutions from the steepest descent needs the prior 
inversion of the system of equations (l44l in order to obtain (/Xj, v{) as functions of the original kinematical 
variables (£i,yi). The correlators are functions of the variables yi) and the virtuality of the jet 
Q = EQq. After setting yi = Y — li with fixed Y = m(Q/Qo) i n the arguments of the solutions (T5TT > 
and (1551) the dependence can be reduced to the following: C% 23 (4,4,4,^) and C^ M (*i , 4, 4 , F). 

3.1 Predictions for the limiting spectrum A « 

In this subsection we give predictions within the limiting spectrum A < 0.5 for charged hadrons mostly 
composed by pions and kaons. 

In Figj3l the DLA (|35T >. MLLA hump approximation from subsection [2/7] and MLLA (f5TT > three-particle 
correlators are displayed, as a function of the difference (£\ — £2) = hi^/xi) for two fixed values of 
£ 3 = ln(l/x 3 ) = 4.5, 5.5, fixed sum (£1 + £ 2 ) = | ln(sia;2)| = 10 and finally fixed Y = 7.5 (virtuality 
Q = 450 GeV and Aqcd = 250 MeV), which is realistic for the LHC phenomenology lfi"3l . The values 
£3 = m(l/x3) = 4.5, 5.5 (#3 = 0.011, X3 = 0.004) have been chosen according to the range of the 
energy fraction Xi -C 0.1, where the MLLA scheme can only be applied and in particular, the range 
x < 0.014, where the cumulant correlator is dominant (see subsection 12. 81 ). 

In Figf4j the DLA d35l ), MLLA hump approximation from subsection 12.7 1 and MLLA dFTT ) three-particle 
correlators are displayed, in this case, as a function of the sum (l\ + £2) = \ \n{x\X2)\ for the same 
values of £3 = ln(l/x3) = 4.5, 5.5, for x± = X2 and Y = 7.5. The range 7.0 < | m(xiX2)| < 13.0 has 
been chosen according to the condition x < 0.014 discussed in 12. 8 1 

As expected in both cases, the DLA and MLLA three-particle correlators are larger inside a quark than 
in a gluon jet. Of course, these plots will be the same and the interpretation will apply to all possible per- 
mutations of three particles (123). As observed and written above, the difference between the DLA and 
MLLA results is quite important pointing out that overall corrections in 0{^/a~ s ) are quite large. Indeed, 
the last behavior is not surprising as was already observed on the treatment of multiplicity fluctuations 
of the third kind, where 11321 

(n(n-l)(n-2)) r . , . 

^ ^ = 2.25 [1 - (1.425 - 0.021n/)Va7] , 

(n) G 

(n(n-l)(n-2)) G 

— =4.52 [1 - (2.280 - 0.018n/)V5£] . 



For instance, for one quark jet produced at the Z° peak of the e + e~ annihilation (Q = 45.6 GeV), one 
has a s = 0.134. Replacing this value into the previous formula for a quark jet multiplicity correlator, 
one obtains a variation from 4.52 (DLA) to 0.83 (MLLA). That is one of the reasons for DLA has been 



18 



Y=7.5 ; 1=0.1 ; |Ln(x x )| 



:10 



Y=7.5;X=0.1 ;|Ln(x x )|=10 



4 

3.75 
3.5 
3.25 
3 

2.75 

So" 2.5 
O 

2.25 
2 

1.75 
1.5 
1.25 





DLALn(1/x 3 )=4.5 

DLA Ln(1/x 3 )=5.5 

Hump Ln(1/x )=4.5 " 

Hump Ln(1/x )=5.5 - 

MLLALn(1/x 3 )=4.5 

MLLALn(1/x 3 )=5.5 











. DLALn(1/x 3 )=4.5 
DLA Ln{1/x 3 )=5.5 
. Hump Ln(1/x 3 )=4.5 
Hump Ln{1/x 3 )=5.5 
. MLLA Ln{1/x 3 )=4.5 
MLLA Ln(1/x )=5.5 



Ln(x 2 / Xl ) 




Figure 3: Three-particle correlations inside a gluon jet (left) and a quark jet (right) as a function of 

£ 1 -£ 2 = ln(x 2 /x 1 ) for i x + £ 2 = I Hxix 2 )\ = 10, £ 3 = ln(l/x 3 ) = 4.5, 5.5, fixed Y = 7.5 in the 
limiting spectrum approximation A ~ 0. 



known to provide unreliable predictions which should not be compared with experiments. From Figj3] 
the correlation are observed to be the strongest when particles have the same energy Xi = Xj for fixed 
Xk and to decrease when one partem is much harder the others. Indeed, in this region of the phase 
space two competing effects should be satisfied: on one hand, as a consequence of gluon coherence and 
AO, gluon emission angles should decrease and on the other hand, the convergence of the perturbative 
series k± = XiEQi > Qo should be guaranteed. That is why, as the collinear cut-off parameter Qq 
is reached, gluons are emitted at larger angles and destructive interferences with previous emissions 
occur. Moreover, the observable increases for softer partons with x% decreasing, which is for partons less 
sensitive to the energy balance. In Fig J4] the MLLA correlations increase for softer partons, then flatten 
and decrease as a consequence of soft gluon coherence, reproducing for three-particle correlations, the 
hump-backed shape of the one-particle distribution. Because of the limitation of the phase space, one 
has < 1 for harder partons. Finally, in FigJ5] we display the three- particle correlators as function 
of the sum | Ixl(xiX2Xs)\, for x\ = x 2 = x 3 ; when compared with Fig]4]and FigfJJ the correlators are 
shown to be larger. That is why, and as expected, the correlations are the strongest for particles having 
the same energy-momentum x\ = x 2 = x%. In these figures, the MLLA hump approximation is seen 
to become larger than the DLA correlator for smaller values of x than those close to the hump region, 
which is unphysical. This is due to the fact that this approximation should not be trusted beyond the 
hump region | £ - Y/2 |< a oc F 3 / 2 , 3Y/2 = 11.25 in this case. 

The MLLA hump approximation from subsection 12-71 is observed to be larger than the MLLA solution 
from the steepest descent of the one-particle distribution but one should bear in mind that this is only 
an approximation made for the sake of clarity in the interpretation of the solutions. In particular, from 
Figj3]one can observe a smoother descent for the slope of the correlators in this case than that given from 
the more exact steepest descent. This difference comes from the role played by the iterative corrections 
displayed in FigJU which decrease the correlators away from the hump region when one of the partons 
becomes harder than the others. Near the maximum Xi = Xj of the correlators, the difference between 
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Figure 4: Three-particle correlations inside a gluon jet (left) and a quark jet (right) as a function of 

h + ^2 = | hi(rciX2)| for x\ = X2, is = ln(l/x3) = 4.5, 5.5, fixed Y = 7.5 in the limiting spectrum 
approximation A « 0. 



the two approaches is O ^^7oj and should decrease for Xi — > 1, according to (l68l i. 
3.2 Predictions beyond the limiting spectrum A ^ 

The approximated evaluation of the one-particle distribution from the steepest descent method made 
possible the evaluation of the two-particle correlations beyond the limiting spectrum approximation, that 
is for Qq / A.qcd- Accordingly, it makes also possible the evaluation of the three-particle correlators 
C r (^1,^2,^3)^) and C n (h,^2,^3,Y) beyond this limit A / 0. This parameter, also known as 
hadronization parameter, guarantees in particular the convergence of the perturbative approach a 8 <l. 
In Fig J6] and Fig|7]we display the same set of curves beyond the limiting spectrum (A = 1.5) as in 
Figj3]and Fig|4]in the limiting spectrum (A ~ 0), with the exception of curves coming from the hump 
approximation. The value of A in this case was evaluated for Qq ~ 1 GeV, which corresponds to the 
proton mass, and Aqc d = 250 MeV. As observed the correlation increases with A and the range where 
C^ 3 ) > 1 becomes larger in this case. 



4 Conclusions 



In this paper we provide the first full pQCD treatment of three -particle correlations in partem showers 
and a further refined test of the LPHD within the limiting spectrum approximation and beyond. The 
evolution equations satisfied by this differential observable have been obtained for the first time and 
the differential version of the equations has been solved iteratively. It has been possible to interpret 
the analytical solution in terms of Feynman diagrams describing the process and to evaluate it from the 
steepest descent method applied to the single inclusive distribution. The correlations have been displayed 
in the range x < 0.014, where the process is dominated by three particles emitted from the same partonic 
cascade following the QCD AO described in FigfJJ and Figf2]i Furthermore, four-particle correlations 
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Y=7.5 ; X=0.1 ; x^x =x Y=7.5 ; ; x =x 2 =x 




3|Ln{x t )| 3|Ln(x 1 )| 

Figure 5: Three-particle correlations inside a gluon jet (left) and a quark jet (right) as a function of 
h + ?2 + ^3 = | In (#i £2 £3) | for x\ = X2 = X3, fixed Y = 7.5 in the limiting spectrum approximation 
A ps 0. 



Y=7.5 ;X=1.5 ; |Ln(x x )|=10 Y=7.5 ;>.= !. 5 ; |Ln(x x )|=10 




Lnfx^) Lnfx/^) 

Figure 6: Three-particle correlations inside a gluon jet (left) and a quark jet (right) as a function of 

h-l 2 = ln(x 2 /xi) for i x +l 2 = | ln(xix 2 )| = 10, £ 3 = ln(l/x 3 ) = 4.5, 5.5, fixed Y = 7.5 in the 
limiting spectrum approximation A = 1.5. 
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Figure 7: Three-particle correlations inside a gluon jet (left) and a quark jet (right) as a function of 

h + ^2 = | hi(rciX2)| for x\ = X2, is = ln(l/x3) = 4.5, 5.5, fixed Y = 7.5 in the limiting spectrum 
approximation A = 1.5. 



have been computed at DLA so as to show that the inclusion of higher order corrections for more than 
three particles would rather be a cumbersome task. The correlations have been shown to be strongest 
for the softest hadrons having the same energy x\ = X2 = X3 in both quark and gluon jets, increasing 
as a function of ln(xj/xj) and | ln(xiXj)\ when Xk softens, that is for partons being less sensitive to the 
energy balance. 

Coherence effects appear when one or two of the partons involved in the process is harder than the 
others, thus reproducing for this observable the hump-backed shape of the one particle distribution. 
Away from the maximum at Xi = Xi, because of limitation of the phase space, one has < 1. 
Predictions beyond the limiting spectrum for heavier charged hadrons as compared with pions and kaons 
show that the correlations should increase as the parameter Qo equals the mass of such hadrons and the 
range where > 1 has been enlarged beyond this limit. The last statement is not surprising because 
soft gluon emission gets suppressed between the two scales Qo an d ^-qcd for A / 0, thus decreasing 
the particle yield inside the whole jet. This measurement would in particular provide an additional 
and independent check of the LPHD for massive charged hadrons. As was shown in 12.41 the DLA 
solution of the evolution equations provide general features of the observable showing its unreliability 
to be compared with the experiment. That is why, the MLLA shape and overall normalization of this 
observable should be compared with the data. In the case of pp collisions at the Tevatron, since diet 
events consist of both gluon and quark jets, in order to compare data to theory, a parameter f g for mixed 
samples of quark and gluon jets was chosen iTTTl . In pp collisions at the LHC, the same procedure can be 
applied so as to measure the two- and three-particle correlations. Furthermore, MLLA corrections have 
been shown to be larger for three than for two particles, that is to increase as the number of particles 
increases. 

As was the case for two particles, the three-particle correlations are larger inside a quark than in a gluon 
jet. Same trends have been observed in HERA and LEP data for soft multi -particle fluctuations in 11341351 . 

Finally, we give the first analytical predictions for intra-jet three-particle correlations in view of forth- 



22 



coming measurements by ATLAS, CMS and ALICE at the LHC. 
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A MLLA approximation 

In (I23ab . for ln(l — z) <C In x and In z <C In x, we perform the following Taylor expansions: 

Q(3) (1 _ z) _ q(S) w ln(1 _ z) *Q® + C( ) (71) 

dii 



Qfhl -z)- QW ) (G k (z) -Q k ) + [ GfMz) - Q\V ) (Q fc (l 



j(2) 



,(2), 



,(2) 



Qfc) 



ln(l - z) 



dQ jf < Gi oa + (c {2) o {2 A dQk 

dh [G k - Q k ) + [G tJ - Q l3 ) ^ 



+ 0(a s 



(Qi - Gi(z)) (Qj(l -z)- QA Q k ps ln(l - z)(Qi - G { 



dQj_ 



Q k + 0(a s ). 



(72) 



(73) 



Since none of these terms contribute to MLLA O(yfal), they will be dropped hereafter. In equation 
d23b| ), we perform the following approximations in the hard fragmentation region, 



G^iz) - G%>) {G k {l -z)- G k ) ps hxz\n{\ - z] 



,(2) 

'a 



jgj? G k 
dh dh 



+ 0(a s ), 



(G l -G i (z))(G j (l-z)-G j )G k 



lnzln(l - z)— -— L G k + 0(a s 



Neither d74l nor (|75T ) contribute to MLLA. The other terms in (123bb can be written as, 

dQ® 



2Q^\z) - G (3) ps (2Q^ - g(3)) + 2 In. 



+ 0(a s 



(74) 
(75) 

(76) 



2 (Qg } (z) - Gf ) (Q k {\ - z) - G k ) « 2 (Qg } - Gf) (Q fc - G fc ) (77) 
+21n(l - z) ( g) - ^ + 21n,(Q fc - G fc )^ + 0(a s ), 

(2Qi(z)Qj(z) - GiGj)G k ps - GiGj)G k + In* f Qi^ + ^T^') + ( 78 ) 

(d - 2Q i (z)){2Q j {\ -z)- Gj)G k ps - 2Q i )(2Q i - G,)G k - 2(2Qj - Gj)G k hxz^G k 
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+ 2(G t - 2Q t ) ln(l - z)^G k + 0(a 



(79) 



such that only the first terms in ([76), (|77T ). (|78T ) and (1791 will be kept in the following. Furthermore, we 
make use of the identity |fT9l 



dz(l - z)$fO) (g@)(z) + 

such that G^{z) - zG^ can be replaced by, 

G {n) {z) - zGW ->(!-«) [G (n) (z) + (G (n) (z) - G (n) )l (1 - z) 



G (n) (z) +ln 



(n = 1, 2, 3) in the r.h.s. of equations (I21bb . (I22bl ) and (I23bl ). Indeed, terms oc In z, ln(l — z) provide 
NMLLA corrections C(a s ) which improve energy conservation; however, their inclusion goes beyond 
the scope of the present paper. 



A.l One and two particle distributions at small x 



The MLLA integro-differential version of equations (121al2 lbb and (|22bl22ab is obtained after integrating 
over the regular part of the splitting functions, such that 



Qi,y = ^ I ' di'^' + y)G i (i',y)-*^(£ + y)G i (£,y), 



Cf 
N eJl) 



4iV c 



G 



hV 



o'rtil' + y)Gi{e, y) - ai l{l + y)G i (£, y), 



(80) 
(81) 



with y%(£ + y) = Me + y+x) , and the two-particle correlations {Af) = A ( ff - A t Aj) EMU, 



G?) 



7 2 (£ + yj )G® (£, yj , Vij ) - ai l{£i + Vj )G?) (4, M,) 



+ (a - b)~f (£i + yj)G(£i,yj + r]ij)G(£i + 



(82) 



(83) 



with ^(ii + yj) 



(I80I8TT) and (I82l83l) with respect to "f , one has El 



, after accounting for hard corrections 0{yfa~ s ). After differentiating 



Q 



N c ,u N c 4 
Gify = 7oG, - a7o(Gj^ - A)7oGi) 

from where the following useful relations hold in MLLA [19], 
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+^(7o 2 : 
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(84) 
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(86) 
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Corrections oc /3 in d84l and (f85T >. which are NMLLA, account for the running of the coupling constant 
a s and those oc |, a, (a — b) account for energy conservation in the hard parton splitting region. The 



MLLA gluon inclusive spectrum is given by the solution of (1851) [ 6 ] and can be written in the form 11141 : 

T(B) n dr 



Gi(£,y) 



P 



o Jo 



-Ba 



7T 



Fb{t,v,£), 



(89) 



1 y 

where the integration is performed with respect to r defined by a = - In — + ir and with 

2 Z. 
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■ sinh a 



+ y a 



Po sinh a 
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I B (2^Z(T,y,£)), 



Z(r,y,i 



+ y 



a 



cosh a 



y-t 



sinh a 



Po sinh a \ y + I 

B = a/(3o and Ig is the modified Bessel function of the first kind. The formula in (|89l corresponds 
indeed to the so-called hump-backed plateau, which describes the energy spectrum of soft hadrons in the 
limiting spectrum approximation Qq = Aqcd |6]|28]. This result is well known and constitutes one of 
the strikest predictions of pQCD. The corresponding solution of (l84l for Qi(£, y) can be obtained from 
(I87T ) with accuracy 0{yfa~ s ). The system of differential evolution equations for two-particle correlations 
follows from ([82]> and d83), such that Q9) 



(2) 
ij 

(2) 
ij 



QiQj 



GiGj 



ly N ( 



— 7o G y 



(2) 



3Cp 
4 iVi 
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-To (G^ - Po-fo 



(90) 
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+ (a - 6) 7o 2 [(GiG,), - Po^GiGj] . (91) 



In Ifl9l . the system ( 1821831 ) was solved iteratively after replacing G 

,(2) 



( 2 > = cg^G, and Qg) 



GQ ijQiQj in dSB and (|90]> respectively. The MLLA solutions of d90]> and (|9B, which are to be used in 
the present paper read [ 19 1 



C 



: (2) 
(2) 



1 - <5f - b (Vv + ^jY«) 
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(92) 



(93) 



which were evaluated by the steepest descent method over the single inclusive distribution in ||25l . We 
have introduced the following notations and functions [ 19 ], 



X 



XI 



InC 
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l 3 



-o(iy, 

d X lJ 
di 



0(7, 



o 2 ), 
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dx 



0y 



O(7o 2 ); 



7o 



Xl (i>i,y + i>j,y) + Xy + i>i,i) = C(7o), 



(94) 
(95) 
(96) 



where, following from (|38l l and (|39l , we have evaluated the corresponding order of magnitude of these 

quantities in powers of the anomalous dimension 70 oc ^/a~ s . The solution is iterative with respect to 

■ (2) 

corrections \ and 0\, which need the prior evaluation of the DLA solution C G . of the equations. 
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B Iterative solution of the evolution equations 



Let us first solve the equation (1501 . For the sake of simplicity, it is much easier to solve the equivalent 
equation: 



&2 =7o 2 ^-a 7o 2 (Gf -/3 7o 2 ^) + (a-6) 7o 2 { 
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One has to substitute the following in the l.h.s. of the equation ([97 
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Thus, after normalizing by 7 2 GiG 2 G3, one finds, 
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while for the other terms in the r.h.s. of the same equation one finds, 
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The r.h.s. provides the following contribution 
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After adding (|98l and (|99l and equating with (1 100 1 ) together with some algebra in between, one finds the 
solution written in (IFTT ). Following the same iterative procedure 

Q (3) = CqLq^Qz, q! 2) = 4 2 i^> G(3) = 4 3 L G i G 2G 3 , 
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for the quark jet evolution equation written in (1491 , one has, 
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Finally by adding and subtracting (e% + e%) in every term oc (Cq . — lj in the l.h.s. of (11011 ) one finds 
455). 



C Steepest descent evaluation: reminder from [25 J 

The evaluation of the integral representation by the steepest descent method at small x -C 1 (or large 
£ 3> 1) and very high energy 7»1 leads to the result, 
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2 -y/7rcoshiADeL4(/x, i/) 



(/i — i/) cosh fi cosh f + cosh \i sinh f — sinh fi sinh f 



sinh fi cosh f 

The logarithmic derivatives of the spectrum given in (1651 and (l66l ) were derived from (1102 1 and it was also 
shown that (11021 ) reproduces the Gaussian shape of the inclusive distribution near the hump £ max ps Y/2. 
From (1 1 02b . one has indeed, 



G(£,y) 



nv%[(£ + y + \) 3/2 - A 3 / 2 ] 
where the MLLA £ max reads, 



1/2 



exp 
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(103) 



Setting a = and A = in the previous expressions one recovers the DLA results, which are needed for 
subsection 12.41 The functions entering as a function of (//, v) in d65l ) and d66l ) are the following, 

cosh sinh /i cosh v — (// — v) cosh z^ — sinh v 



Q{n,v) 
KM 

CM 



(fi — v) cosh \i cosh v + cosh /x sinh v — sinh /u cosh v ' 

1 {a — v) cosh /j — sinh /j 

Slllll v 

2 (/j — v) cosh ^ cosh v + cosh /X sinh v — sinh /x cosh v ' 
3 1 (/x — v) cosh z^ sinh /x + sinh z> sinh fi 

2 ^ 2 (/x — i>) cosh /i cosh f + cosh /x sinh v — sinh /x cosh v ' 



L(/x, z/) + tanh v coth /i (1 + K (/j, z/)) . 
The expressions for the two particle correlations follow from (|92l ) and d93l [25], 

,(2) l-6 70 (ew+e«)-^ 



C. 



1 + 



1 + 2cosh(/Uj - fj,j) + A'(m,Vi,iJ,j,i>j) + Sf 



C, 



(2) 



n — - 



CV — 1 H — (b — a)7o — 

G v 2 V ; ,u l + cosh(/x 



(104) 

(105) 

(106) 
(107) 

(108) 
(109) 
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where, 
and 



$i = Polo 



2 sinh 2 [ &-£L 



3 + 4smh 2 (^ 



Q(m,Vi) + Q(iij,vj)), 



(110) 



-070 e Mi + — sinh(//j — [ij)(Qi — Qj) + cosh fix tanh v 2 + cosh fi 2 tanh z^i 

— sinh ^ tanh z/» coth — sinh fij tanh coth ^ 
+ sinh(/ij — fij) ^ tanh i/j coth /UjQi — tanh Z7 coth fijQj 

— Polo cosh/ij — sinh/ijCj + cosh^- — sinh^Cj + sinh(^j — fJ-j)(CiQi — CjQj) 
+ cosh /ij tanh ^'(1 + JCj) + cosh /ij tanh ^(1 + Kj)] . (HI) 



The solutions (11081 ) and d 109b are the ones to be used in this paper for the evaluations of the three-particle 
correlations and will be directly inserted in the solutions (TSTI ) and (|55T > respectively. 



C.l Corrections , and ei, ei 



For the computation of these corrections, one only needs to take the DLA part of the logarithmic deriva- 
tives of the one-particle distribution tpn = 7oe Mi and ip^ y = 7oe _/il , such that after replacement in (|54ct 
and d58c| ) one finds, 



If 



1 

To 
1 

7o L 



X?' (e^ 1 + + e^ 3 ) + (e m + e^ 2 + e 
xf (e" w + + e^ 3 ) + xf (e^ + + e 



where 



X* = A)7oy 



2 tanh 



x/ 



tanh^ 



+ 2 cosh(^, — /ij) 2 



A 2 ) 



with 



C. 



(2) 
Gy 



1 + 



+ 2cosh(/ij - /ij) 
1 



A 2 ) 



(2) 



c 

1 + 2 cosh (//i - ' y y 



1 + 



Ci? 1 + 2cosh(/ij — fjLj) 



Accordingly, replacing ip^ i = 7oe Mi and ipi >y = 7oe Mi in (I54el) and (158el) . one has 

ei = 1 [0 (e"^ 1 + e"^ 2 + e^ 3 ) + ( y (e w + e^ 2 + e" 3 )] , 
7o 



1 

ei = — 

To 



C £ (e"^ 1 + e"^ 2 + e"^ 3 ) + 4 (e^ 1 + e^ 2 + e 



i-i-.i > 



(112) 
(113) 

(114) 
(115) 
(116) 

(117) 
(118) 



where Q, Q and ( y , C, y should be found from the DLA expression of written in (135T ). for Ca = N c in 
a gluon jet and Ca = Cf in a quark jet. Introducing the parametrization in (/u, 1/), one has respectively, 



j(3) 
'G123 



|(2) 
'G12 



5(2) 
'G13 



j(2) 

'G23 



(119) 



1 

+ r 



2 1 + cosh(/ii - fj, 2 ) + cosh(/ii - ^3) + cosh(// 2 - ^3) 
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1 

+ - 



2 1 + cosh(/ii - fj, 2 ) + cosh(/ii - ^ 3 ) + cosh(^ 2 - m) ' 



and 



4(3) 
''Q123 



1+ c, 



i(2) 
'Q12 



1 + c. 



;(2) 



1 + C. 



i(2) 
'Q23 



(120) 



+ 



+ 



4(2) 
-O12 



1 + c, 



a(2) 



1 + c, 



4(2) 



1 



2Cj? 1 + cosh(^x - /i 2 ) + cosh(^! - /x 3 ) + cosh(/i 2 - ^3) 

Nl 1 

2Cp 1 + cosh(^i - m) + cosh(^i - fi 3 ) + cosh(/i 2 - ^3) 



Thus, in order to get Q and one should start from (II 1911201 1 and make use of 

dfM _ dm 
dy dy 



dm d/ij 2 e^Qi - e^Qj dm dfij 2 e ^Qi - e 

■¥"if--wo o ' -^-^r-PWo o 



Therefore, everything is ready for the computation of 







1 /i(3) 



C. 



1 



c. 



(3) 



(3) ^Gra/' s f ^(3) °Gi23,!/' 







G123 



c, 



1 /5(3) 

(3) °Ql23, 



C. 



1 /5(3) 
4(3) °Qi23,3/ - 



(121) 



Q123 



For instance, 



4(3) 



12/4(2) 13/7W , *J/i 

X| + ^ L G!3 + Xt L 



13/4(2) 



53/5(2) 



G13 



G23 + 2 1 + cosh(/xi - /x 2 ) + cosh(/ii - /x 3 ) + cosh(/x 2 — // 3 ) 



y 12/5(2) +v 13/5( 2 ) + v 23/j(2) 
Xt L Gl2 + ^ L Gi 3 + Xt L-G23 



(2) 
G12 



1 + c. 



5(2) 
'G13 



1 + c. 



4(2) 
'G23 



2 [1 + cosh(/xi - /z 2 ) + cosh( / ui - /x 3 ) + cosh(/j 2 - // 3 )]' 



sinh( / ui - (i 2 ) 



+ sinh(/ii - /i 3 ) 
1 



9£ d£ 



+ sinh( / u 2 - /x 3 ) 



_ dm 
di dl 



1 



2 [1 + cosh(/xi - /x 2 ) + cosh( / ui - /x 3 ) + cosh(/j 2 - m)Y 



sinh( / ui - /i 2 ) 



dm _ dji2 
di dl 



dm _ fy±2 

di dl 



+ sinh(/ii - m) 



dm 9 m 



di di 



and 



,5(3) 



(>W _ ^12/5(2) -13/5(2) , -23/5(2) , gig --c gig - -c V2 3 

°Qi23,^ « ^ Ql2 -r Xt ^ Ql3 f ^q 23 ^ 2C7 F 1 + cosh^i - /i 2 ) + cosh^ - /x 3 ) + cosh(/i 2 - m) 



+ sinh( / u 2 - m) 



dm _ dm 

di dl 



(122) 



xfc^xfc^ z+ xfq 



4(2) 
?23 



iV, 



4(2) 



1 + c, 



4(2) 



1 + c, 



•4(2) 



2Cf [1 + cosh(/ii - m) + cosh(/ii - m) + cosh(/i 2 - |U 3 )] Z 



+ sinh(/xi - m) 



dm dm 



dl 



+ sinh(/i 2 - m) 
1 



dm dm 



dl di 



2C|, [1 + cosh(m - m) + cosh(/ii - m) + cosh(/i 2 - m)Y 



sinh( / ui - /i 2 ) 



sinh(^i - m) 



dm _ dm_ 

dl dl 



dm _ dm_ 

dl di 



+ sinh(^i 



M3j 



dm dm 



+ sinh(/i 2 - m) 



dm dm 



dl dl 



(123) 



For derivatives with respect to y, it is enough to replace i by y in the previous expressions. In FigJU we 
display €\{l\, £ 2 , -£3, Y) as a function of the sum | ln(xix 2 )| and the difference {l\ — i 2 ) = ln(x 2 /xi) 
for two fixed values of £ 3 = ln(l/x 3 ) = 4.5, 5.5, x\ = x 2 and fixed sum (l\ + i 2 ) = \ ln(xix 2 )| = 10 
and fixed Y = 7.5. As expected, this correction decreases the correlations away from the hum region 
and for harder particles. 
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C.2 Hump approximation 



In this approximation, we consider the energy of the three partons to be close to the maximum of the 
single inclusive distribution | t — Y/2 a oc Y 3 / 2 for i = 1, 2, 3. In ('25], it was demonstrated that, 



k~Y/2 ^ 1 2l , ti,j~Y/2 



j (i + m + -n { ), ip, 



1 41 



7o(l - fa + -fx]) 



3 y 



(124) 



2y + f 

for a,/?o,A = 0, which is DLA. In the same approximation one has the following for a,/3o / and 

A = 0, 

ti,j~Y/2 

Aij ' « 2-\-{ni-Hj) - 070(2 + fii + fj,j) - 2/3 7o, 

where 



0* - Mi) 



2 * J , 



Y 



ii,j~Y/2 



3 1 



y 



(125) 
(126) 



Moreover, 



V ~ o/3o7o(^i - Hj) = Wolo 



Y 



(127) 



since 70 I - ± T J - 



y 



, 8\ was neglected in this approximation. 



Applying the previous expansions to (153aH53dl ) and (!57aH57dl ). it is easy to find: 



AT, 



(2) 



Qij 



0, 



N, 



(3) 



3c 



G 



5 _ | ln(xix 2 x 3 )| \ l_ 



Po \2 HQ/Qo) J VHQ/Qo) 



' " 3C| 2 Y 



D§ ] =J Dg } =8+9 



3a 
Wo 

8+9 



5-2 



4-4 
Y 



m(x 2 /a7) 

HQ/Qo) 

| ln(xix 2 x 3 )|\ 



ln(x 3 /xi) 

HQ/Qo) 
i 



+ 9 



ln(x 3 /x 2 ) 
HQ/Qo) 



6/3 



y/&>HQIQo) 



(128a) 
70, (128b) 

,(128c) 



HQ/Qo) J ^HQ/QoY 



+9 



y 



-6/3 7o-3a 5-2 



?i+l 2 +4T 

y 



1 7o, 



30 



(2) 



(t. 



C 



(2) 



c, 



(2) 



D 



(3) 
Q 



N. 



(3) 



36 /5 | ln(xiX2X3)| \ 



ft V2 ln(Q/Q ) J ^HQ/Qo) 



I, ( g _ g | ln(xja:j)| 



1 + 



ln(Q/Qo) y y/ln(Q/Q ) 



1 



3 + 9 



ln(gi/gj) 
MQ/Qo) 



1 — 6(5 — 3£i±£i 



1 + 



HQ/Qo) Vfo 
— i To 



5-3 



HQ/Qo) J y/HQ/Qo) 



3 + 9(^y 




C (2) . 






C (2) _ 

Gij 





2/3 7o-a 5-3^p 70 



9 + 9 



ln(x 2 /xi) 



1 + 4 (? 

2 

+ 9 



lxi X^Xj I 



VHQ/Qo) 



70, (128d) 
(128e) 

(128f) 
(128g) 



a) 70 (5-3 £4+ ^ 



y 



H x 3/xi) 



HQ/Qo 
9a ( 5 | ln(xiX2X3)| \ 



ln(Q/Q ) 
1 



+ 9 



ln(x 3 /x 2 ) 

HQ/Qo) 



6/3 



y/PoHQIQa) 



(128h) 



/3 V2 ln(Q/Q ) / VMQTQo)' 

2 / « « \ 2 



V Y 

(72 U Gl23 



+9 



+9 



Y J \ Y 
3a /5 | ]xv{x\X2X^)\ \ 



-6/3 7o-9al J 70, (128i) 



1 



/ft V2 ln(Q/Q ) 7 y/HQIQo) 
3a /5 ^i + ^ + M 1 



(128j) 



Vft V2 ln(Q/Q ) 7 y/HQ/Qo) 

D DLA solution of the 4-particle correlations 



Below, we display the expressions related to subsection |2.9i In the l.h.s. of the evolution equation d69l ), 
we define 



i(4) 
■■1234 



,(4) 
1 1234 



\^ 3 -A 1 A 2 A 3 )A i -(A {3) 



134 



A1A3AAA2 



1 (3) 
1 234 



A 2 A 3 AAAi 



(129) 



1 (3) 
l 124 

i(2) 
•■14 " 



- I AW - A\A± ) ( A 



(2) 
23 



A\t - A X A 2 A± ) A 3 - ( 4? - MA 2 ) ( Afl - A 3 A A 



A 2 A 3 ) + (Af 2 ) 



1 (a) 

i 34 



1(2) 

L 13 



AiA 3 )A 9 A 



12/14 



1(2) 



A1A4 )A 2 A 3 + A 



1 (2) 
*23 



A 2 A 3 ) >M 4 + (4? " ^2^4) Ax A 3 
+ - A3A4J A X A 2 - AiA 2 A 3 A±. 

In the DLA solution ( PTOl ) of the equation d69l ), we have introduced the expressions: 



Hi 



C 



(2) 
12 



1 + c 



Ho 



C 



(3) 
123 



+ (4? 



4(2) 
-13 



+ (C (3) 

r I u 124 



1 + c. 



4(2) 
'23 



1 + c 



;(a) 

'34 



1 + c 



4(3) 
'134 



1 c 



(2) 



12 



c 



(2) 
13 



Die 



1) + (c. 

(2) 
21 



;(3) 

'234 



1 + c 



4(2) 
-14 



2 C 



4(2) 
'12 



-2 fd 



(2) 
23 

(2) 
13 



2 C 



4(2) 
-11 



2 C 



4(2) 
'23 



2 C 



4(2) 
'21 



2 C 



4(2) 
'34 



(130) 
(131) 
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H S = 1 + (cl - l) + (cg» - l) + (eg> - l) + (eg - l) + (eg' - 1) (eg» - 1) (B2) 
+ (4? - 1) (eg' - 1) + (eg? - 1) (4? - 1) - (4? - 1) - (cg» - l) 

- (4? - 1) - (4? - 1) - (4? - 1) - (4? - 1) ■ 
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